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Abstract- Data from two different satellites, a digital land cover map, and digital census 
data were analyzed and combined in a geographic information system to study the effect of 
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urbanization on photosynthetic productivity in the United States. Results show that 
urbanization can have a measurable but variable impact on the primary productivity of the 
land surface. Annual productivity can be reduced by as much as 20 days in some areas, 
but in resource limited regions, photosynthetic production can be enhanced by human 
activity. Overall, urban development reduces the productivity of the land surface and those 
areas with the highest productivity are directly in the path of urban sprawl. 
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I. INTRODUCTION 

As we begin to recognize the scope of human influence on Earth’s ecosystems, it is important to 
understand how specific forms of human induced land transformation affect the dynamics of 
Earth's biological systems. Land transformation due to human activity has taken many forms 
historically starting with fire management, herding practices, the development of agriculture, and 
culminating with urbanization including industrial development [1-3]. Past studies of human 
impacts to the biosphere estimate that between one third to one half of the planet's land surface 
has been transformed by human action [4] and that nearly 40% percent of the yearly products of 
photosynthesis are appropriated by human beings [5]. Recently, more attention is being paid to 
urbanization as it constitutes one of the more ecologically disturbing land transformation 
processes and urban areas are expanding rapidly as human populations grow in size, affluence, 
and technological capability [6-9]. 

While there are numerous issues that could be addressed concerning urbanization, in this study 
we focus specifically on how urbanization affects the seasonal dynamics of photosynthesis and 
total yearly primary production. We examine transformed land cover types around large 
metropolitan areas, compare them to their lesser-urbanized or non-urbanized counterparts, and 
we make a continental scale assessment of productivity for three classes of urbanization. The 
relative impact to photosynthetic production is measured by observing the change in the 
normalized difference vegetation index (NDVI) as a function of the degree of urbanization. 
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II. METHODS 

We used a combination of satellite data, information extraction techniques, and geo-spatial data 
from map and tabular sources to carryout the analysis. Nighttime images from the Defense 
Meteorological Satellite's Operational Linescan system were used to create a thematic map 
portraying the extent and spatial distribution of urbanized, peri-urban, and non-urban areas in the 
coterminous United States. A digital land cover map generated by the United States Geological 
Survey (USGS) was used to define “pre-urban” land cover types and photosynthetic productivity 
was estimated using monthly maximum normalized difference vegetation index (NDVI) values 
from AVHRR data. The data sets were co-registered and staged in a geographic information 
system (GIS) to intersect and analyze the parameters in a spatially explicit way (Fig. 1). Using 
the DMSP/OLS urban thematic map as a template, NDVI values were calculated within each 
land cover type both inside and outside urbanized areas. The difference in primary productivity 
between the urbanized and non-urbanized portions of each land cover type was calculated by 
integrating the NDVI values over time. Land cover specific analyses were carried out by 
examining the land surface parameters within 200 x 200 km plots (40,000 km 2 ) for seven major 
cities in the coterminous United States. The approximate centroid of each city was located in the 
center of the sample plot so that urban, peri-urban, and non-urban areas were included in each 
plot around the cities. 

A continental scale view of the impact of urbanization on photosynthesis was achieved by 
grouping the land cover classes together and comparing the different NDVI based productivity 
measurements for urbanized, peri-urban, and non-urban lands for all of the pixels in the 
coterminous U.S. 
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A. Mapping Urbanized Lands from Space 

There is no internationally standardized definition of "urbanized land" and map data representing 
this form of land cover at regional or super-regional scales are rare. As such, there is need of a 
methodology for identifying urban land use and mapping it synoptically. We used nighttime 
images from Defense Meteorological Satellite Program's Operational Linescan System to create 
the urban map for this study. Originally designed to map moon-lit cloud cover for nighttime 
aircraft navigation and weather forecasting for the United States Air Force, the DMSP/OLS 
operates at extreme sensitivity collecting image data at a moderate spatial resolution (2.7 km 
pixel) across a broad vis-NIR band (.4 to 1.1 pm). The resulting images show a dramatic picture 
of urbanization through the detection of city lights during cloud free acquisitions [10]. ' 

Previously, DMSP/OLS data have been used to estimate population [11] and indicate energy 
consumption [12,13]. 

The DMSP/OLS data used in this study were collected from October 1, 1994 - March 31, 1995, 
between 20:30 and 21:30 local time. The data were screened for cloud cover and ephemeral light 
sources, and re-projected onto a 1 km grid to conform with other global data bases [14]. These 
are composite data containing dozens of "hits" for each grid cell over the six month long 
acquisition period. The DN for each 1 km cell is a percentage representing the frequency of 
illumination for that cell (i.e. the number of times the cell was registered as lit divided by the 
sum of the OLS acquisitions for that cell). Using previously developed thresholding techniques, 
reasonably accurate maps showing urban core areas, peri-urban, and non-urban lands were 
generated based on the percentage of time a grid location was registered as lit during the data 
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collection period (See Imhoff et al. [15] for a detailed description of this process and accuracy 
assessment). 

A thematic product identifying three classes of urban land use was generated from the 
DMSP/OLS data: 1) Urbanized, 2) Peri-urban, and 3) Non-urban land. Urbanized areas 
exhibited high levels of nighttime light emission (>89%); Peri-urban areas exhibited 
considerable but unstable illumination (5-88%); and Non-urban lands were never observed as 
being lit (<5%, 5% being considered noise [14] indicating no or little development. These classes 
and their areas compare well with census data estimating urban and non-urban areas in the 
United States [15]. Demographic parameters such as population and housing density were 
indexed to the DMSP/OLS map product by overlaying it with digital census data in the GIS 
[16,17]. Using the aggregate average (i.e., all seven urban study plots combined) the Urbanized 
land class had an area- weighted average population density of 1064 persons and 445 housing 
units per square kilometer, the Peri-urban Land class had 100 persons and 38 housing units per 
square kilometer, and the Non-urban Land class had an average of 14 persons and 5 housing 
units per square kilometer. A table showing the average population and housing densities 
associated with each of the urban areas featured in this article can be found in Table I. 

B. Pre-Urban Land Cover Determination 

For part of our assessment of the impact of urban land transformation on photosynthetic 
production we constrained the analysis by comparing changes within a land cover class. This was 
done to limit the possibility that the land use change signal would be lost or confounded by 
differences between land cover classes. We used the United States Geological Survey's North 
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America Seasonal Land Cover Regions land cover classification to define land cover types for 
our analysis. This map is part of the North American Land Cover Characteristics Data Base and 
cover types are identified, in part, by their temporal spectral characteristics using 1-km AVHRR 
data. The Seasonal Land Cover Regions Map was chosen because it defines 205 land cover 
classes for the United States providing greater resolution than it's IGBP or Global Ecosystems 
counterparts [18-20]. However, the classification is hierarchical in nature and all of the Seasonal 
Land Cover types can be identified in the IGBP system. The North America Seasonal Land 
Cover Regions land cover classification is well suited for this analysis because it defines a large 
number of classes based primarily on vegetation type but contains no distinct "urban" land cover 
class. The large number of land cover classes and the absence of a distinct urban class provided 
the best possible "pre-urban" land cover view for our analysis. 

C. Estimating Photosynthetic Production via Satellite 

Data from the National Oceanic and Atmospheric Administration's polar-orbiting advanced very 
high resolution radiometer (AVHRR) sensor can be used to produce a vegetation index directly 
related to chlorophyll density contained on land [21-26]. The normalized difference vegetation 
index (NDVI) is expressed as [(Channel 2 - Channel 1) / (Channel 2 + Channel 1)], where 
Channel 1 records reflected light from 0.55-0.70 pm and Channel 2 records reflected light from 
0.72-1.1 pm. Previous studies have shown that the index is an accurate estimate of absorbed 
photosynthetically active radiation [27,28] and is strongly correlated to biomass [29-31]. 

In this study, a 1 km AVHRR composite data set of monthly maximum NDVI acquired from 
April 1992 through March 1993 was used to estimate photosynthetic production. The 


6 

11 


Imhoff et al. 2000 


compositing process minimizes atmospheric effects, scan angle effects, cloud contamination, and 
solar zenith angle effects [32]. By integrating over time, monthly maximum NDVI data from 
AVHRR composites have been shown to be useful for estimating the annual and semi-annual 
primary productivity of vegetation on land [33-36]. 

To assess the impact of urban land transformation on productivity, a three-way merge was 
accomplished in the GIS where the land cover map was overlain on the NDVI composite, and 
the DMSP/OLS urban map was applied as a final layer. Every pixel in the merged data set had 
fifteen attribute values associated with it: 

1) urbanization class from the DMSP data, 

2) land cover type from the USGS land cover map, and 

3) maximum NDVI values for each of the twelve months of the composite year from the 
AVHRR data set. 

A mean urbanized and non-urbanized NDVI value was then computed for each land cover class 
(and for each month) by averaging all of the pixels inside and then outside the urbanized areas 
defined by the DMSP/OLS layer. The average monthly maximum NDVI values were multiplied 
by the number days in the month to produce units of NDVI*days. Summed over the year this 
measure has been used to represent annual primary production (see Prince [34] and Prince and 
Goward [35] for a more detailed description of this method of evaluation). A temperature 
threshold, using the thermal channel on the AVHRR (Channel 5, 1 1.5-12.5 pm), was used to 
identify pixels where the mean monthly temperature was at or below 0° C. NDVI's for these 
pixels were assigned a value of 0 as relatively little photosynthesis occurs at these temperatures. 
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A set of NDVI*day signatures covering a period of one year were then generated for each cover 
type showing the differences inside and outside the urbanized zones. 


III. RESULTS 

A. Urban Impact Signatures 

By combining the three data sets in a GIS we were able to generate comparative signatures of 
photosynthetic productivity for each land cover class as they exist inside and outside the 
urbanized areas as defined by the nighttime DMSP/OLS data. The results indicate that the impact 
of urbanization on photosynthetic production depends on two primary factors: the type or density 
of development, and the climatic/physiographic environment. Although the analysis was 
performed on over 100 cover types existing in and around seven large metropolitan areas in the 
US, four primary patterns of the seasonal dynamics of the impact of urbanization on productivity 
appeared repeatedly in the data (Table II). 

1) Continuously lower photosynthetic production in the urbanized portion of the land cover type 
(Fig. 2), 

2) Extended growing season in the urbanized areas resulting in higher productivity at the 

beginning and end of the growing season and lower peak season productivity (Figs. 3 and 4). 

3) Higher productivity in the urbanized areas throughout the year (Fig. 5), and 

4) Higher productivity in the urbanized areas during a resource limited part of the seasonal cycle 
but lower yearly summed productivity (Fig. 6). 
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Continuously lowered productivity in urbanized areas was found for numerous cover types 
around all the metropolitan areas researched in this study. This effect is common for heavily 
urbanized land cover types and is well represented by a Savanna land cover type found around 
Miami, Florida (Fig. 2). This land cover type has been heavily used for urban development in 
south Florida. With the exception of deep winter, the productivity of the urbanized portion of this 
land cover class was reduced more or less uniformly throughout the year resulting in a summed 
loss of 22.4 NDVPdays for the year. This pattern is most likely the result of intense urbanization 
where a substantial portion of the vegetation in the urbanized area has been converted to 
impervious surfaces supporting roads and other infrastructure. These surfaces are nearly devoid 
of plant species and will substantially lower the monthly maximum NDVT calculated for these 
areas [37]. 

The second pattern, which is characterized by an extended growing season in urbanized areas, 
was found to be common for numerous cover types in seasonally cool or cold parts of the United 
States. In fact this pattern was the most commonly repeated signature for all of the cover types 
we studied. Urbanized areas within these cover types show an earlier start of photosynthesis in 
the spring, a lower peak production, and a slower reduction at the end of the growing season. 
Occasionally mid-winter photosynthetic activity was also observed in these areas possibly due to 
a heat island effect. The yearly impact of urbanization in these areas is variable, some 
demonstrating a net loss, others showing a net gain in productivity. 

Examples of this pattern can be seen in two Mixed Forest cover types around Chicago, Illinois. 
The forest cover type in Fig. 3 (Mixed Forest #1) is composed largely of Aspen, Maple, Oak, 
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Pine, and Spruce and the forest cover type in Fig. 4 (Mixed Forest #2) is composed of Maple, 
Beech, Birch, and Pine. In both cases, productivity in the urbanized portion of the land cover 
class is enhanced from March through June, reduced relative to the non-urbanized areas from 
July through September, and slightly higher again in the winter months. Both exhibit rapid 
declines in productivity with the onset of Autumn (October-November). We interpret these 
trends to result from a combination of urban heat island effect and the cultivation of grasses and 
shrubs in lower-density urbanized areas. This compensated for the negative effect of 
urbanization on productivity if development is not too intense. Higher temperatures would 
prolong the growing season for the grasses and shrubs in this cover type while the reduced peak 
season is again due to higher percentage of impervious non-vegetated surfaces. 

However, the net yearly effect of urbanization on the productivity of the two forests is quite 
different: the Aspen-Maple-Oak forest showed a yearly loss of productivity (5.56 NDVI*days) 
while the Maple-Beech forest showed an increase in productivity (4.29 NDVI*days). The rapid 
fall off of productivity for both urban and non-urban zones in autumn indicated a substantial 
population of deciduous trees and shrubs in both areas. We speculate the higher urban 
productivity in the winter could indicate a composition change in vegetation cover type, from 
mixed forest to a higher incidence of planted grasses and evergreen shrubs. 

The third and fourth patterns are similar to one another in that both occurred in semi-arid, arid, or 
seasonally-arid regions where the urbanization is not too intense. These patterns show an 
enhancement of productivity inside the urbanized areas, especially during the warmer months, 
and sometimes throughout the year. We attribute the increased productivity to irrigation and the 
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alteration of natural land cover by cultivation of introduced plant species. For example, the 
Open Shrubland land cover type in and around Denver, Colorado exhibited higher productivity 
in the urbanized area as compared to natural land cover through most of the year indicating 
conversion of natural semi-arid land cover to one of irrigated grasses, shrubs, and trees (Fig. 5). 
In this case, the urbanized portions of this land cover class gained a yearly sum of 10.73 
NDVI*days over the non-urbanized areas. 

A fourth "Mediterranean" pattern was found for land cover types around urban centers along the 
West Coast of the United States which experience cool-wet winters and hot-dry summers. 

Higher productivity was found in the urbanized parts of the land cover class during the dry 
summer season but during the wet winter season the productivity in the non-urban areas greatly 
exceeds that of the urbanized part. This pattern is well expressed in a Woody Savanna land cover 
type around Sacramento, California (Fig. 6). During the months from May-December, the dryer 
warmer season, the urbanized areas in this land cover class gained 4.86 NDVI*days over the 
non-urbanized area. However, in the cool wet season, from January through April, the non- 
urbanized portions of the land cover type gained 8.9 NDVI*days over the urbanized areas. The 
yearly summation showed a loss of 4.0 NDVI*days due to urbanization. The increased 
productivity of the urbanized areas in the dry season is thought to be a result of irrigation. 
However, the rapid fall off and overall reduced productivity of the urbanized area in the winter 
wet season could be due to various factors. An increased proportion of impervious surfaces in the 
urbanized areas could explain this pattern but the replacement of native evergreen species with 
exotic deciduous species of trees and shrubs, which go dormant in the winter (wet) season, could 
also contribute to this pattern. 
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B. Continental Scale View 

We created a continental scale view of the impact of urbanization on photosynthesis by 
calculating the average monthly maximum NDVI for all the land area in the coterminous United 
States for the three classes of urban land as defined by the DMSP/OLS data. The analysis 
resulted in a quantification of the seasonal pattern of productivity for all of the vegetation land 
cover types in Urbanized, Peri-urban, and Non-urban areas (Fig. 7, Table III). This aggregate 
approach revealed the relationship between land use in the United States and photosynthetic 
production. The land surface with the highest productivity was found in the Peri-urban areas. 
The summed yearly production for the Peri-urban land area was 135.3 NDVPdays, or 10 
NDVI*days longer than the Urbanized Lands (125.3 NDVI*days for Urbanized Lands). The , 
Peri-urban lands in the U.S. contain some of the most productive soils [9] and are frequently still 
in agricultural use or remain as secondary forest. These areas support the highest photosynthetic 
productivity in the United States and are directly in the path of urban expansion. 

Perhaps the greatest surprise was the low productivity calculated for lands in the Non-urban 
class. The Non-urban class had a yearly summed productivity of 1 17 NDVI*days - that is 18.3 
NDVPdays shorter than the Peri-urban class and 8.3 NDVI*days shorter than the Urban class. 
This is primarily because the Non-urban class, defined by a nearly complete absence of 
illumination at night, contains vast areas of resource-limited lands in deserts, high altitude 
montane environments, wetlands, beaches, and rocky or degraded land. Human activity is 
limited in these areas by the same resource limitations that drive other biological systems. It is 
interesting to note that, despite the aggregation of land cover types in the continental scale 
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analysis, the time behavior of photosynthesis for the Urbanized land and Non-urban land was 
very similar to the extended growing season dynamic seen for various individual cover types as 
described above. 


IV. CONCLUSIONS 

We quantified the effects of urban land transformation upon photosynthetic production through a 
multi-level, geo-spatial analysis of nighttime and daytime satellite data, land cover maps, and 
ancillary census information. The effects of urban land transformation on productivity are 
variable, depend on the degree of urban land transformation and the climatic/physiographic 
settings. In general, urbanization appears to lower the overall photosynthetic productivity of the 
landscape. While irrigation and urban heat islands may offset the effects of increasing the area of 
impervious surfaces, the overall trend is toward reduced yearly photosynthetic productivity in the 
urbanized zones. Viewed at the continental scale, urbanization resembles other biological 
processes: growth originates where resources are concentrated; uses available resources; and 
diminishes where resources are limited. The clear trend is that the most productive land (peri- 
urban land) is directly in the path of urbanization. On average, we can assume that the 
transformation from peri-urban to urban will result in an annual loss of about 10 NDVI*days of 
photosynthetic production. Overall, about 3% of the land surface of the United States is now 
urbanized and about 15% of the land surface is in the peri-urban class. What the overall affect of 
the continuing process of urbanization will be on the aggregated photosynthetic production 
potential of the continent remains to be considered. 
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TABLE I 

MEAN POPULATION AND HOUSING DENSITY FOR URBAN RELATED CATEGORIES 

DERIVED FROM THE DMSP/OLS DATA 



Population Density (persons/km 2 ) 

Housing Density (housing units/km z ) 

Miami 



Non-urban 

3.5 

1.4 

Urbanized 

1034.6 

477.2 

Chicago 

Non-urban 

18.9 

7.0 

Urbanized 

1263.4 

490.1 

Denver 



Non-urban 

4.7 

2.2 

Urbanized 

923.9 

441.9 

Sacramento 



Non-urban 

13.2 

5.7 

Urbanized 

1147.8 

461.7 

USA Total * 



Non-urban 

13.7 

5.4 

Peri-urban 

100.0 

38.3 

Urbanized 

1063.7 

445.1 


* USA Total is the aggregated average of 200 x 200 km plots around seven cities (Denver, CO; 
Chicago, IL Miami, FL; Sacramento, CA, San Francisco, CA; Seattle, WA; Washington D.C.). 
Demographic values from 1990 census. 
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TABLE II 

TEMPORAL SUMS OF NDVI (NDVPDAYS) FOR FIVE LAND COVER TYPES AROUND MAJOR 
URBAN AREAS IN THE UNITED STATES. 


Pattern 

IGBP Land Cover 

1992 









1993 






Apr 

May 

Jun 

Jul 

Aug 

Sept 

Oct 

Nov 

Dec 

Jan 

Feb 

Mar 

SUM 

I 

Savanna 

Miami, Florida 
Urbanized 

11.76 

11.04 

9.30 

13.27 

13.36 

12.79 

12.21 

11.31 

13.02 

11.47 

11.45 

13.24 

144.22 


Non-Urban 

12.51 

14.35 

13.50 

15.10 

15.72 

15.42 

15.13 

13.92 

14.88 

11.16 

10.39 

14.54 

166.61 

II 

Mixed Forest #1 

Chicago, Illinois 
Urbanized 

7.99 

16.27 

17.62 

19.31 

17.98 

18.38 

17.67 

4.94 

6.45 

6.87 

5.14 

5.89 

144.52 


Non-Urban 

7.32 

13.27 

16.21 

20.25 

19.72 

20.13 

17.45 

4.55 

5.00 

6.13 

5.11 

5.11 

140.22 

II 

Mixed Forest #2 
Chicago, Illinois 
Urbanized 

9.05 

15.65 

15.99 

17.85 

17.16 

17.26 

17.01 

5.19 

6.82 

6.13 

5.06 

6.36 

133.16 


Non-Urban 

8.57 

13.71 

14.90 

20.29 

20.29 

20.50 

17.22 

5.25 

5.80 

5.89 

6.16 

6.51 

138.57 

in 

Open Shrubland 
Denver, Colorado 
Urbanized 

6.34 

6.71 

10.01 

9.05 

8.93 

8.81 

6.94 

6.26 

4.01 

5.13 

4.15 

4.19 

80.52 


Non-Urban 

5.35 

5.36 

9.02 

7.96 

7.96 

7.95 

8.05 

5.31 

2.16 

3.61 

2.95 

4.11 

69.79 

IV 

Woody Savanna 
Sacramento, California 
Urbanized 

13.38 

9.77 

8.07 

8.43 

7.42 

7.31 

7.09 

7.65 

4.37 

9.72 

7.96 

12.01 

103.17 


Non-Urban 

17.09 

9.93 

7.47 

7.31 

6.21 

6.08 

5.94 

7.30 

4.99 

10.98 

9.71 

14.19 

107.20 


Statistics were gathered within 200 x 200 km plot around each city. 

Mean monthly maximum NDVI values were calculated inside and outside the urbanized areas 
for each cover type. A cold mask based on AVHRR channel 5 data was used to factor out pixels 
with temperature Odegress C or less. 


TABLE III. 

TEMPORAL SUMS OF NDVI (NDVI*DAYS) FOR ALL LAND IN THE UNITED STATES DIVIDED 
INTO THREE CLASSES OF URBANIZATION BASED ON DMSP/OLS DATA. 



1992 









1993 




DMSP Classification 

Apr 

May 

Jun 

Jul 

Aug 

Sept 

Oct 

Nov 

Dec 

Jan 

Feb 

Mar 

Sum 

Urbanized Land 

10.16 

12.38 

12.51 

13.10 

12.61 

12.46 

11.99 

8.83 

8.64 

8.05 

6.85 

7.76 

125.35 

Peri-urban Land 

10.41 

13.24 

14.09 

15.35 

14.65 

14.23 

12.67 

8.86 

8.78 

8.53 

6.90 

7.59 

135.29 

Non-urban Land 

9.16 

11.43 

12.45 

13.53 

12.26 

11.91 

10.95 

7.53 

7.68 

7.92 

5.80 

6.41 

117.03 


Mean monthly maximum NDVI values were calculated for each urban class. A cold mask based 
on AVHRR channel5 data was used to factor out pixels with temperature 0° C or less. The 
urbanized, peri-urban, and non-urban lands represents approximately 3, 15, and 82% of the total 
surface area of the continental U.S. respectively. 
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Imhoff et al. 2000 



Figure 1. The geospatial data used in this study. (Top) An urban classification map based on 
DMSP/OLS satellite data [15] showing three classes of land use: urbanized land (red), peri-urban 
land (yellow), and non-urban land (black). (Middle) The United States Geological Survey’s 
North America Seasonal Land Cover Regions land cover map [ 18-20]. (Bottom) An NDVI 
image from the AVHRR satellite monthly maximum NDVI data set [32,33]. A single month 
(July) is shown as an example. The urban map, the land cover map, and the AVHRR data set 
were co-registered for analysis in a GIS. 
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Savanna - Miami, Florida 



Month 


Fig. 2. Profile of photosynthesis for urban and non-urban areas in and adjacent to Miami, 
Florida. Paved surfaces result in lower aggregated photosynthetic production in the urbanized 
parts of the land cover class. 
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Mixed Forest #\ - Chicago, Illinois 



Fig. 3. Productivity signatures for urbanized and non-urbanized areas of a Mixed Forest cover 
type in the vicinity of Chicago, Illinois. This pattern was found most typically in seasonally cold 
geographic regions. Photosynthetic production starts earlier in the urbanized areas, has lower 
peak growing season productivity, and has higher production in the winter months. An urban 
heat island effect and the ecosystem alteration are probably responsible for the increased 
productivity in the urbanized areas which counter balances the negative effects of impervious 
surfaces. In this example, however, there was still a negative overall impact on productivity due 
to urbanization. The urbanized part of this land cover class had a net loss of 5.56 NDVI*days 
relative to the non-urban part. 
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* Mixed Forest #2 - Chicago, Illinois 



Month 


Fig. 4. Seasonal pattern for another Mixed Forest cover type around Chicago, Illinois. The 
pattern is very similar to one shown in Fig. 4, however, in this case, the yearly total productivity 
is greater in the urbanized area (net gain of 4.29 NDVI*days). The increase is attributed to the 
urban heat island and a higher incidence of planted grasses and evergreen shrubs in these areas. 
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Open Shrubland - Denver, Colorado 



Month 


Fig. 5. Open Shrub Lands around Denver, Colorado show overall enhanced photosynthetic 
production in urbanized areas. The greater productivity in the urbanized portions of this cover 
type is probably due to replacement of native cover type with species supported with irrigation. 
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Fig. 6. Woody Savanna cover type around Sacramento, California showing a Mediterranean 
seasonal pattern. Irrigation and ecosystem alterations with exotic species may both contribute to 
the urban signal. Note higher productivity for urbanized areas in the peak dry months. This is 
almost certainly due to irrigation. In the winter- wet season, the exotic deciduous species in the 
urbanized areas go dormant reducing the NDVI of these areas. 
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Continental United States 
All Cover Types 



Fig. 7. Annual NDVI based productivity signatures for the three DMSP derived classes of 
urbanization for all of the surface area of the continental U.S. The value for each month is the 
mean NDVI*days calculated for all land cover types in each class of urbanization. 
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